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An analysis of the transient heat transfer between a thin circular tube 





and the incompressible fluid moving through the tube is made for the case 
where the temperature of the inlet fluid is kept constant.^^Bpth radial con- 
duction of heat in the wall and zhe heat loss at the outer surface of the 
cylinder are taken into consideration. It is shown to be especially easy to 
caJLculate temperature of both fluid and tube in the initial period by means 
of the attached figures. 

So far it is common to use the temperature of the wall calculated by the 
assuniption that it is constant radially. It is shown that for the insulated 
tube in the initial period such temperature is equal to the temperature of 
the outer surface of the tube. The temperature of the inner surface of the 
tube may be appreciably different from such temperature even for the metallic 
tube. The difference is extremely large for the tube made from insulating 


material . 
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F(a_,7;z) 

G = pV 

Gi(C^t)^G2(^,t) 

Iy(n) 

k 

ko 

ki 

K 

Kv(t]) 
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Li(t]) 

Ni(t]) 

X(&) 

£-1(5) 

Ku 

P 
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Pr 


specific heat at constant pressure of inner fluid, 
kcal/kg deg C 

heat capacity per unit volume of tube material, kcad/m^ deg C 
thickness of the wall 

dimensionless functions defined by eq.uations (24) 
confluent hyper-geometric function 
dimensionless functions defined by equations (45a) 
dimensionless functions 
mass flow per unit area, kg/m^ hr 

non-dimensional functions defined by equations (45b) 
modified Bessel function of the first kind and of vth order 
thermal conductivity of tube material, kcal/m hr deg C 
thermal conductivity of fluid kcal/m hr deg C 
dimensionless constant, h/RUo 
dimensionless constant, k/PUi 
dimensionless constant, pcpk/2CUoR 

modified Bessel function of the second kind and of vth order 
length of the tube, m 
dimensionless constants 
dimensionless function 
dimensionless function 
Laplace transform of 6 
inverse Laplace transform of 5 
Nusselt number, 2RQo/kf, dimensionless 
variable of Laplace transform 
pressure of the inner fluid, atm 
Prandtl number, t^Cp/kf, dimensionless 



3 


R 

Re 

t 

Ti(x,t) 

Ta(x,r,t) 

Ti 

Ts 

Uo 


Ux 


V 

X 

y 

a 


a* 


A 

3 

3* 

B 

7 

5 

6 


radial distance, m 
inner radius of the tube, m 
Reynolds number, 2 RG/h, dimensionless 
time, hr 

temperature of inner fluid, deg C 
temperature of tube, deg C 
inlet temperature of inner fluid, deg C 
temperature of outer fluid, constauit, deg C 

film coefficient of heat transfer between inner fluid and tube 
surface, kcal/m^ hr deg C 

film coefficient of heat transfer between outer fluid and tube 
surface, kcal/m^ hr deg C 
velocity of inner fluid, m/hr 
axial distance, m 

S 

dimensionless constant, (l + a)Ui/UQ 
dimensionless constant, (m/K) - 1 
dimensionless constant, a(l + a) 
dimensionless constant, a* + 1 
dimensionless constant, A + 1 
dimensionless function defined by equation (l4) 

(is - Tg)/(Tj_ - Ts), dimensionless 
X(6) 


t L|, dimensionless 

2LUc,l/RGCp, dimensionless 
T] i’/R, dimensionless 
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For any intermittent or blowdown hypersonic wind tunnel it is required to 
establish an air flow at elevated values of stagnation temperature. Thus^ air 
is preheated in a heater and led to the test section through a settling chamber 
and nozzle. For the design of such components of the wind tunnel it is 
required to determine the transient temperature of the wall of each component 
as well as the temperature of air. This problem may be reduced to determine 
the transient heat transfer between a circular tube and the incompressible 
fluid flowing through the tube, since the main flow field is subsonic and the 
effect of compressibility can be neglected. 
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Such a problem has been treated already by several investigators. Rizika[l]^ 
calculated the transient temperature of fluid taking into account the heat loss 
at the outer surface of the tube. His result is inaccurate aind does not contain 
the analysis of the teniperature history of the wall. Judd[2] calculated tran- 
sient temperature of both fluid and wall for the insulated tube. Vfe extended 
Judd’s ceQ-culation so that the effect of heat loss at the outer surface of the 
tube was included[3]« In these treatments it was assumed that the temperature 
of the wall is constant radially. The use of such an assumption may be justi- 
fied if the thermal conductivity of the material of tube is large. Indeed, 
this was verified experimentally by Judd for a copper tube. 

For the high temperature and/or high pressure application such as a hyper- 
sonic wind tunnel, however, we are forced to use such material as stainless 
steel or ceramics. For these materials thermal conductivity is smaller and 
the validity of the assumption is doubtful. Therefore such an assumption is 
replaced by a more plausible assumption of the thin wall in this paper. 

ANALYSIS 

The following assumptions are imposed on the solution: 

1. The temperature ' of inner fluid is function of time and axial distance 
from the inlet only. 

2. The temperature of tube is function of time, axiel distance from the 
inlet, and radial distance from the axis. 

3" The temperature of outer fluid is constant. 

4. The effect of thermal conduction is negligible in the axial direction. 

5- The velocity of inner fluid is axial, constant, and uniform. 

6. Both the inner and the outer radii of the tube are constant, and the 

ratio of the difference between these radii to the inner radius is small. 

^Numbers in brackets designate references at end of paper. 
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7 • l''Iaterial constants do not depend upon temperature . 

8. There are no energy sources within the tube material itself. 

9 - The film coefficients of heat transfer between the fluid sind tube are 
uniform and constant over the inner and outer tube surfaces for a constant 
fluid mass flow rate. 

10. The effect of radiation is negligible. 

11. The ratio of the thermal capacity per unit volume of air to that of 
tube material is negligible. 

Consider the system shown in Fig. 1 . Let an incompressible fluid be flow- 
ing in a circular pipe with a flow rate G; the temperature of this fluid at 
the axial distance x from the inlet and at time t is Ti(x,t). The tempera- 
ture of the circular pipe is T2(x,r,t)^ where r is the radial distaince from 
the axis. The temperature of the fluid outside the pipe is assumed constant 
and equal to Tg. 

Equating the sum of the thermal energy crossing the boundary of an ele- 
mental section of fluid gives 

3ll , , ,, 2U(J , 

o=p -57 + p=p -§i- - = - [Ti - TsU.K.t)! 


Similarly^ for an elemental section of the wall, one obtains 

p dr ^^2 k 3t 

Now we can neglect the third terms on the left-hand sides of equations (l) 
and (2) by assumption 4 . Adopting the nondimens ional variables, these equa- 


^ + K II = 6UA,tp) - 0 ( 3 ) 

8^6 1 dB dB 


tions are reduced to 
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where 


K = £££ k 
^ ~ C 2 UoR 


Since k/(2UQR) = 0(l), K is negligible by assumption 11. Therefore 
equation ( 3 ) reduces to 

|| + 0 = s(^,i,q>) 

Equations (4) and (5) are fundamental equations for this analysis. 
Initial conditions are given by 

e(|,0) = &(|,Ti,0) = 0 


(5) 


Boundary conditions are 


0(0, cp) = 1 


^ *^11=1 


V 


e - 5(i,i,cp) 

6(1,1 + a,cp) 


T)=i+a 


After the following Laplace transformations 

_ poo 

0(1, P) = / 0(|,<p)e”^dcp 

» p 00 

6(l,T),p) = / &(|,T),cp)e'^acp 

*^o 

and using equation (6), equation (4) is reduced to 

8n2 1 01 

The general solution of this equation is given by 

B(l^n>o-) = gi(|,a)lQ(aTj) + g2(l>o,)Ko(aTi) 


(6) 


(7) 


(8) 
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v/here a = and gi; and g 2 are dimensionless functions to be determined 
later. Using this equation, equation (5) is reduced to 

• + 0 = + Ko(°-)S2(l^a-) (9) 

Boundary conditions, equations (j), are transforined to 


0(0, p) = p‘^. 


= 0 - 6(1,1, a) 

\ i/r)=i 

-ki(§^) = S(IA + a,a) 

V T\ = 1 -)-n. 


T) = 1-HX 


Substituting equation (8) into equation (lOc) gives 

^ Ko(3) - kiaKi(p) 

■ Io(13) + kialii 3 ) 22(1.°') 

Using equations (8) and (ll), equation (lOb) reduces to 


(iOa) 

(10b) 

(10c) 


( 11 ) 




rTCo( 3 ) -kioKi(3) ^ K^(P)- kiaKi(3) 

|lol3)+ kialiipT ■ Iol3) + kiaIi(3) lo(a) .+ ..Ko(°-) 


Putting these expressions for gi and g 2 into equation (9) gives 

^0 7- 

- 70, 


(12) 


( 13 ) 


where 
. 7 = 


1 


L Io(a)K^(3) - Io(3)Ko(0') - kia)l^(a)Ki(3) + Ii(3)Kp(a)} 


-1 


( 11 ) 


^ " k^a Ii(a)Ko(3) + Io(3)KiCa) - kia|Ii(a)Kxi3) - Ii(.3)Ki(,a)| J 

The solution of equation (13), satisfying the boundary condition (lOa), is 

0 ( 1 , P) = ( 15 ) 


given by 
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Substituting equations (il) and (l2) into equation (8) gives 


5(1, 11, a) = [Ko(3)Io(c^Ti) - Io(3)Ko(a7i) - kia<Ki(3)Io(aii) + Ii(3)Ko(aTi){ ]e/y 

(16) 


where 

t = Ko(p)lo(a) - Io(3)Ko(a) - kia jKi(3)Io(a) + Ii(3)Ko(a)f 

-koa[Ix(a)Ko(3) + Io(3)Ki(a) - kiajli(a)Ki(3) - Ii(3)Ki(a)( ] (17) 


In general it is not easy to get the inverse Laplace transformations of 
equations (l5) and (16). Therefore we restrict our considerations to the cases 
where the wall of the tube is thin and the ratio of the thickness of the wall 
to inner radius of the tube is far less them unity: a « 1. Since 3= (l+a)a, 

one has 

I„(8) = (1 ^ ar (l)Vn(-) 

n=o 

00 Yi 

K,{p) = (1 + a)” (-1)" 

n=o 


Therefore, neglecting the terms of 0(a®), one obtains 

Io(a)Ko(3) - Io(3)Ko(o-) = Ii(a)Ki(3) - Ii(3)Ki(a) = . | (2 - a) 
Io(a)Ki(3) + Ii(3)Ko(a) = (l - a) ^ 

Ii(ot.)K^(3) + Io(3)Ki(a) = i + I a2 

Ii(a)Ki(3) - Ii(3)Ki(a) = . | (2 - a) 

Substituting these expressions into equations (l4) gives 




> (18) 




V+^J 


7 = L 1- 


(19) 
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Vfhere 


L = 


2ki + (l - kx)a 
2ki + |Ui/ko) + 1 - kija 


< 1 


M 


2kp+2ki+2(l- ki)a> (l - 2ki)ag 
koa[2ki + |Ui/ko) + 1 - ki JaJ 


N = M - 


|2ki + (1 - ki)a|a 

Thus equation ( 15 ) is rewritten as 


e(l^p) = e 


-L| 




Now one has 


’ll 


expi 


p + M 


= e 




whose derivation is given in appendix. Therefore one obtains.^ 

e(|,cp) = e'^.^|e-^Io (2 >yi^) + M 

Combining the equations (16) and (18) together with the following 

Ko( 3 )lo(o<il) - Io(p)Ko(a7]) = - | {2a - a^ - 2(r\ - l) + (tj - l)^j 
Ki( 3 )lo(o-n) + Ii(p)Ko(°'l) = (1 - a + a^) ^ + I fa - (n - l)p 


where the terms of O(a^) are also neglected gives 


M. 


where 


fi(n) = 


ja - (n - 1)1‘ 


(ko/ki)a2 + ko(2a - aS) + a^ 


faCl) = 1 - 


2ko + 2{r\ - 1) - (ti - 1)‘ 


How one has 


2ko + 2ki + 2(1 - ki)a - (l - 2 ki)a 2 
^ [p(p + M) ®^^Cp+ m)] ^ ® ^Io(2 >/^)dcp 


- fx(ri) 


( 20 ) 

( 21 ) 

( 22 ) 

equations 

( 23 ) 

(24a) 

(24b) 

(25) 
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whose derivation is also given in appendix. Therefore the inverse Laplace 
transformation of equation (23) is obtained as 



6 (l^n^ 9 ) = f“i(T])e(l/ 9 ) + Mf2(Ti)e-L| / 6 -^Iq( 2 LN|cp)d9 (26) 

*^0 

Using ^ = L|, T = Wcp, A = (m/N) - 1 , B = A + 1 = M/N, equations ( 22 ) and (26) 
are reduced to 

n T 

(27) 




ana 


= tx{Ti)e-^Mi,r) + B ^fi(n) + f2(n)|j e'^MC.T)dT (28) 


respectively, where 


= e■^■'^Io(2^/TT) 


Since fi(l) = 1 - L, fi(l + a) = 0 and Bf2(l) = L, one obtains 
6(C,1,t) = (1 - L)e-^Mi,r) + (B + L - BL) J 


where 


bU,l + a,r) = Bf2(l + a) / ( C,T)dr 

o 


f2(l + a) = 


2ki(l - a + a2) 


2ko + 2a - a2 + 2ki(l - a + a^) 


(29) 

(30) 

(31) 

(32) 


Equations (3O) and (3I) give the temperature of inner and outer surface of the 
tube wall, respectively. Moreover one gets 

5 ( 0 ,ti,t) = fi(Ti) + f2(Tl)(l - e-S-^) (33) 

which gives the temperature distribution in the tube wall at entrance section. 

DISCUSSION 


First it will be shown that the solutions obtained in the previous section 
tend to the steady solutions when t tends to infinity. Here the steady solu- 
tions are the solutions of the equations ( 4 ) and (5) where 86/89 is assximed to 


vanish; 
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0 s(t) = 


(3^) 

(35) 


where 


K = 


ko 


^ i ~ ^a + 


Neglecting the terms of O(a^) and considering 0 <ri-l<a«l give 


K| = 


2kol 


2ko + 2ki + 2(1 - ki)a - (l - 2ki)a2 


■0 


and 


1-K- — lnT) = l-K| 


1 + 




|(n . 1) . I (11 - i) 2 j- 


= fx(il) + fa(il) 


Thus one gets 


and 


®s = exp |- (^1 - C 


&s = {fi(n) + faCn)} exp|- (^1 - 0^j- 


(36) 


(37) 


Now tending t to infinity in equation (27) gives 

^ r°° n 

e(^,oo) = Be-^ / e‘®%(2>/T^)dT 


Since 


J e‘®’^Io(2N/TT)dT = ye"^ ^Io(y)<^ 


where y = 2 \j and b® = B/4^, and 

^"ye-»Vi^(y)ay . ^ axp(^) 

[see reference 3, equation (35) this equation is reduced to 

9(t,eo) = exp|- (1 - 5) ?]■ 


(38) 
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In a similar manner, from equation ( 28 ), one obtains 

S(^,T).co) = |fi(Ti) + f2(Ti)J exp|- (l - ( 39 ) 

Comparing the equations (36) and (37) with equations (38) and ( 39 ) > i't is clear 
that the present solutions tend to the steady solutions when t tends to 
infinity. 

Next it will be shown that the present solutions tend to the solutions 
obtained by assuming that k = oo[3] when k tends to infinity. When k 
tends to infinity one obtains 


L = 1 


B = 


1 - a + a^ ^o 
A a* 


Ui , Ui 

+ 1 - (1 + a) + 1 

Ur 


= p* 


(1 - a + a2)2Uot 


2Uot 


~ a (2 - a)CR ~ a (2 + a)CR 


( 40 ) 


fi(n) = 0 

Bfsin) = 1 

where * denotes the symbol used in reference 3 * Thus it is proved that the 
present solutions coincide with the previous results [equations (25) and (26) 
of reference 3]^ that is 

rn* 

e = + p* / e-®'*’l*y*dTi* 


pT]* 

5 = / 


where = e"^"'H*Io(2-s/Tii*)^ "to the order of a. 
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INITIAL SOLUTIONS 


Solutions for At « 1 will be called as "initial solutions" in this paper. 
Neglecting the te-rms of 0{k^r^) , equations (27), (28), (3^)^ (3I) are 

reduced to 

dii,r) = Fi(S,t) - AF2(C,t) - A^G2ii,r) (4l) 

5(C>n>'r) = {fi(T)Fi(C,T) + f2(il)Gi(C,T)f 

- A[fi(T])F2(^,T) - f2(T))jGi'(C,T) 

- Ga(^T)}] - A^{fi(T)) + f2(n)}Ga(C,T) (42) 

S(C,1,t) = {(1 - L)Fi(C,t) + LGi(C,r)} 

- A[(1 - L)F2(C,t) + LG2(^,t)] 

- A2(i - L)G2(^,t) ( 43) 



^If T < 1, it is easily shown that A^G2(5 ,t) = O(A^t^), so that terms 
proportional to A^G2 (^^t) can "be eliminated from these equations. 
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Graphs of the functions Fx, F 2 , Gi, and Gg have been published in that 
report and presented here also (as Figs. 2(a) through 2(h)) for the reader's 
convenience. Each term of the solutions can be easily calculated by means of 
these figures. 

If there is no heat loss at the outer surface of the thin tube, one 
obtains Ui = 0 and kx = 00, so that A = 0,^ f 2 (l + a) = 1 and 

L = [1 + a/{ko(2 - a)J ] ^ (46) 

Therefore temperatures of inner fluid and of outer surface of wall of insulated 
thin tube are equal to those of inner fluid and of tube wall, respectively, 
which are calculated by the assumption of radially constant wall temperature. 
Temperature of inner surface of wall of insulated thin tube is given by the _ . 
linear combination of these temperatures: 

6(^,1, r) = (1 - L)Fx(^,t) + LGx(C,t) (4?) 

where L is given by equation (46). 

The first order effect of heat loss at the outer surface of the wall 
might be easily taken into consideration for many cases of practical importance 
For this purpose, equation (42) is rewritten as 

= /fx(l)Fl(C;T) + Bf2(T))Gx(^,T)J 

- Ajfx(ii)F2(C,T) + f2(n)Ga(C>T)} - A^{fi(r]) + f2(n)} Ga( ^^ t) 

Usually the second and the third terms are far smaller than the first term. 

In such cases this equation is reduced to 

&(C^T>t) = fx(il)Fi(C^T) + Bf2(il)Gx(C^T) (48) 

^It follows that each first term of equations (4l) through (44) represents 
the solution for the insulated tube , and the second and the third terms 
represent the correction for the heat loss. 
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In a similar manner one obtains 

(49) 

= (1 - L)Fi(^,t) + LGi(C,t) (50) 

&(^,1 + a,x) = Bfa(l + a)Gi(C,r) (51) 


Thus there is no appreciable effect of heat loss on the temperature of inner 
fluid. If a and a/ki are small quantities, ^■/ko = 0(l) and the second order 
small quantities are neglected, one gets 


Bf2(l + a) = 1 - a/(2ki) 

L = 2 ^ 1 - ki 

2 + a/ko (2 + a/kQ)^ 


(52) 

(53) 


SAMPLE CALCULATION 

* 

As examples temperature distributions of both wall material and air have 
been calculated for a settling chamber of a hypersonic wind tunnel. Here 
three kinds of material are considered as possible tube material for such a 
wind tunnel: stainless steel (l8Cr-8Ni), alumina brick and alumina castable. 

Dimensions of the tube are: 

' length of the tube Z = 3 di 

inner radius of the tube R = O.325 m 

thickness of the wall d = 0.04l m 

so that a = d/R = 0.126 . Uq was calculated from the empirical relation[4] 

Nu = 0.023 Re°‘®Pr°*^ 

for the fully developed turbulent flow in a circiiLar pipe where Nu = 2RUo/kf, 
Re = 2 RG/p, Pr = pcp/kf and p is the viscosity of air. Ui was calculated 
by means of the enpirical formula [4] 
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0.25 


/^T > — T \ 

Ui = 0.929 f — I 4; , kal/m^ hr deg C 


where <Ts> is the mean value of the temperature of the outer surface of tube. 
Since <Ts> was not known a priori, it was determined by a trial and error 
procedure. Constants for inner air are given by Table 1, where M is design 
Mach number and P is pressure. Constants for the settling chamber are 
given by Table 2, where ^/x, r/t, and all dimensionless quantities are the 
mean vailues of four which correspond to the four cases in Table. 1. The vari- 
ation for A, L, and fz{l + a) for these four cases is less than 7 percent. 
Calculated temperatures of inner air and of inner and outer surfaces of tube 
wall are shown in Table 3 both inlet (^ = 0) and exit (C = Ci) except the 
temperature of inner air at inlet which is equal to unity. Since At < 0.021, 
initial solution was used for the calculation. The temperatujse of tube wall 
is highest at inlet. This temperature was also calculated by the former 
method[3] where the temperature of tube wall was assumed constant radially. 

The results are included in Table 3 sis s(0,t) and &(0,t)*). From equation (40) 
one obtains 

^ 2 - a 

■ T)* = T 

(2 + a)(l - a + a^) 

so that the difference 'between t]* and t is O(a^) . It may be clear from 
Table 3 that 6(t) is greater than in the amount of O(a^). Thus it 

could be concluded that the thin-wall approximation estimates the temperature 
of tube wall higher than the accurate value. From the designing point of view 
this overestimate is in safety side. 

The temperature of inner air was also calculated by the former method 
for C = ^1* This agrees with the present result in three figures for all 
cases. This confirms equation (49)* It might be seen that Tables 2 and 3 
also confirm equations (50) through (53)- 
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The temperature distributions in the tube wail are presented in 
Figs* 3 through 5* I't will be seen in the refractory materials temperature 
decreases rapidly near the inner surface and then decreases gradually to the 
temperature of the outer surface. 


CONCLUSIONS 

There is an appreciable temperature gradient in tube wall even if the 
material of the tube is such metal as stainless steel. This temperature 
gradient is very large if the materiel is such insulating refractory materials 
as alujnina brick and alumina castable, in which case teniperature decreases 
rapidly near the inner surface and then decreases gradually to that of outer 
surface . 

Temperature of tube well calculated by the assumption that this temperature 
is constant radially is equal to the temperature of outer surface of the tube 
in the initial period, provided a/ki « 1. It might be dangerous to use 
such temperature for the design since the temperature difference between 
^ inner and outer surfaces is sometimes very large. 

Temperature of inner fluid in the initial period calculated by the thin- 
wall approximation is the same as that calculated by the assumption of 
radially constant wall teniperature. 
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APPENDIX 

1. Derivation of equation (21) 

One has, from reference page 244, 

eVP - 1 = i[ Ii(8nW)] 

and from reference page 129, 

P(p + M) = £ e'^f(cp) (54) 

where F(p) = X[f(cp)]. Combining these equations gives 



By the convolution theorem, one obtains 

^ ^ ^ ■ ^}] " e-^Io(2,/l^) - 1 + M 

Since l/p = X[l] and 



one gets equation(2l) . 

2. Derivation of equation (25) 

One has, from reference 5^ page 245, 

p-ie^P = £[Io(2^/A^)] 
Combining this and equation (54) gives 

p + M ®^p + m) =xj^e”^Io(2^/I^)^ 

By the convolution theorem equation (25) is derived# 



TABLE 1.- CONSTANTS FOR IFl'IER AIR 
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TABLE 2.- CONSTANTS FOR THE SETTLING CHAl-IBER 
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TABLE 3.- TEMPERATURE DISTRIBUTIONS 
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Fig. 1 
Fig. 2(a) 

Fig. 2(b) 
Fig. 2(c) 
Fig. 2(d) 
Fig. 2(e) 
Fig. 2(f) 
Fig. 2(g) 
Fig. 2(h) 
Fig. 3(a) 
Fig. 3(b) 
Fig. 4(a) 
Fig. 4(b) 
Fig. 5(a) 
Fig. 5(b) 


FIGURE CAPTIONS 

Thin circular tube. 

-t 

Temperature functions at ^ = 0: Fi = 1, Fa = 0^ Gi = 1 - e , 

Ga = 1 - (1 + C)e"^. 

Temperature ftinctions at | = 0.01: Fi(O.Ol^O) = 0.99005. 

Temperature functions at | = O.O3: Fi(0.03,0) = 0.97045. 

Temperature functions at | = 0.1: Fi(0.1,0) = 0.90484. 

Temperature functions at | = O.3: Fi (0.3,0) = 0.74082. 

Temperature functions at | = 1: Fi(l,0) = 0. 36788. 

Temperature functions at 1=3: Fi(3,0)= 0.049787’ 

Temperature functions at | = 10: Fi(l0,0) = 0.000045. 

Wall temperature distributions for stainless steel tube at inlet. 
Wall temperature distributions for stainless steel tube at exit. 
Wall temperatixre distributions for alumina brick tube at inlet. 
Wall temperature distributions for alumina brick tube at exit. 

Wall temperature distributions for alumina castable tube at inlet. 
Wall temperature distributions for alumina castable tube at exit. 
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Figvire 3(Q') Nisiki Hayasi and Kenji Inoue 
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Figux’e 5(b) Wisiki Hayasi and Kenji Inoue 


